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Abstract The Diffuse Supernova Neutrino Background (DSNB) is the weak glow of MeV neutri- 
nos and antineutrinos from distant core-collapse supernovae. The DSNB has not been detected 
yet, but the Super-Kamiokande (SK) 2003 upper limit on the i/^ flux is close to predictions, 
now quite precise, based on astrophysical data. If SK is modified with dissolved gadolinium to 
reduce detector backgrounds and increase the energy range for analysis, then it should detect 
the DSNB at a rate of a few events per year, providing a new probe of supernova neutrino 
emission and the cosmic core-collapse rate. If the DSNB is not detected, then new physics 
will be required. Neutrino astronomy, while uniquely powerful, has proven extremely difficult - 
only the Sun and the nearby Supernova 1987A have been detected to date - so the promise of 
detecting new sources soon is exciting indeed. 
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1 Introduction 

The Diffuse Supernova Neutrino Background (DSNB) is the flux of neutrinos and 
antineutrinos emitted by all core-collapse supernovae in the causally-reachable 
universe; it will appear isotropic and time-independent in feasible observations. 
(Hereafter, neutrinos means neutrinos and antineutrinos, and supernovae means 
core-collapse supernovae, unless specified.) It results from the quasi-thermal 
MeV neutrino spectrum emitted by newly-formed neutron stars convolved with 
the rapid redshift evolution of the supernova rate. The cosmic energy density in 
neutrinos from core-collapse supernovae is comparable to that in photons from 
stars, ~ 0.01 eV cm~^, and is ~ 10 times less than that of the cosmic microwave 
background. Some earlier literature referred to "supernova relic neutrinos," but 
"relic" often caused confusion with the 10~^ eV neutrinos from the Big Bang. 

Why is detection of the DSNB an important and relevant experimental goal? 

• Understanding supernovae is crucial to astrophysics and physics. 
Supernovae are central to the cosmic history of stellar birth and death; the pro- 
duction of chemical elements, neutron stars and black holes, cosmic rays, and 
gravitational waves; and to exploiting the extreme physical conditions of core 
collapse to probe the properties of neutrinos and hypothetical particles. New 
observations and syntheses are needed. 

• We cannot understand supernovae without detecting neutrinos. 
The optical supernova reveals little about the collapsing core itself, which is em- 
bedded in the stellar envelope; the total energy in photons is much less than that 
in neutrinos, leaves over a vastly longer time, and depends on the properties of 
the progenitor star. Neutrinos are emitted from the core, revealing its properties. 

• Detecting bursts of neutrinos from nearby supernovae is difficult. 
Prodigious emission and large detectors partially compensate the tiny detection 
cross section. For Milky Way supernovae (distance Z) ~ 10 kpc; 1 pc = 3.1 x 
10^^ cm), a neutrino burst would be detected easily, but the burst rate is only a 
few per century. For supernovae in nearby galaxies (D ^ 1-10 Mpc), where the 
total burst rate is above one per year, much larger detectors will be required. 

• The DSNB is a guaranteed steady source of supernova neutrinos. 
While the probability of detecting even a single neutrino from a distant supernova 
is inflnitesimally small, the number of supernovae per year is astronomically large, 
yielding a nonzero rate. The total rate and spectral shape of the DSNB are new 
probes of supernova neutrino emission and the cosmic core-collapse rate. 

The DSNB has not been detected yet, but the discovery prospects are excellent. 
Super-Kamiokandc (SK) is large enough to detect the DSNB at the level of a 
few events per year. These events are hidden by detector backgrounds, which are 
already low and could be substantially reduced with added gadolinium to detect 
neutrons, a proposed upgrade for which there is active research and development. 

Detection of the DSNB will yield science that is difficult to obtain otherwise. It 
will measure the average supernova neutrino emission spectrum, and comparison 
to burst detections will probe the star-to-star variation. It will also be sensitive 
to the optically-dark failed supernovae expected from high-mass stars. 

This review is written for experimentalists, observers, and theorists in nuclear 
physics, particle physics, and astrophysics, and covers the many aspects of the 
DSNB. It is intended to be an accessible overview of the main themes needed to 
predict, detect, and interpret the DSNB, with details given in references. 
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2 Impossible Dreams of Neutrino Astronomy 

Neutrino astronomy is a unique tool for answering some of the oldest and most 
important questions in nuclear and particle astrophysics (1). What happens 
when a massive star runs out of fuel? Are observed high-energy gamma ray 
sources powered by cosmic-ray protons or electrons? What is the origin of the 
highest-energy cosmic rays? Are there astrophysical neutrino sources unseen 
with electromagnetic observations? Are there new properties of neutrinos, or 
new particles, that can be revealed only with neutrino detectors? 

The immense discovery potential of neutrino astronomy arises from the small 
interaction cross section of neutrinos with matter, which allows neutrinos to es- 
cape from and thus to reveal the interiors of dense and distant astrophysical 
objects, undeflected, undegraded, and unobscured. This powerful magic comes 
at a terrible price, as it is much less likely that the neutrinos will interact with a 
detector than with an entire star. As far as we know, neutrinos interact only via 
the weak nuclear force and gravity (but with only tiny masses). 

In 1934, shortly after Pauli postulated the existence of the neutrino, Bethe 
and Peierls wrote that "If [there are no new forces] - one can conclude that 
there is no practically possible way of observing the neutrino" (2). It took until 
1956 for neutrinos to be detected, by Rcincs and Cowan, using a nuclear reactor 
source (3), for which Reines shared in the 1995 Nobel Prize (Cowan died in 1974). 
This detection allowed a first consideration of neutrino astronomy in principle. 

Prescient remarks about this potential were made long ago, e.g., Refs. (4-8). 
As Bahcall wrote in 1964 about solar neutrinos, before their detection, "Only 
neutrinos, with their extremely small interaction cross sections, can enable us to 
see into the interior of a star..." (9). The detection of neutrinos from the Sun 
and the nearby Supernova 1987A (SN 1987A) shows that the required sensitivity 
has begun to be reached. Davis and Koshiba shared in the 2002 Nobel Prize for 
these first detections in neutrino astronomy, which led a precise understanding 
of the solar fusion reactions, the discovery of neutrino mixing, and confirmation 
that Type II supcrnovae result from the deaths of massive stars. 

No other astrophysical neutrinos have been detected, despite decades of effort. 
In Bahcall's 1989 book. Neutrino Astrophysics, which was almost completely 
about solar neutrinos, he wrote that "The title is more of an expression of hope 
than a description of the book's contents. ...the observational horizon of neutrino 
astrophysics may grow ... perhaps in a time as short as one or two decades" (10). 
Indeed, the first detections of distant neutrino sources now seem imminent, due 
to more certain predictions, based on new astronomical and laboratory data, and 
due to astounding improvements in detector sensitivity. Any detections will be 
of enormous scientific importance, certain to provide new information. 

This review is about the exciting prospects for new results in extrasolar MeV 
neutrino astronomy; as described below, these build on a long history of work in 
experiment, observation, and theory. There are similarly encouraging prospects 
at the TeV (10^ MeV) scale (11,12) and the EeV (10^^ ^gv) scale (12,13). 

3 Framework: DSNB Detection Spectrum at Ecirth 

Here the framework to calculate the DSNB detection spectrum is given, along 
with a sTimmary of its three ingredients: the supernova neutrino emission (Section 
4), the cosmic supernova rate (Section 5), and the detector capabilities (Section 
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6), followed by a return to the framework (Section 7), then SK observations 
(Section 8), background reduction (Section 9), and conclusions (Section 10). 



3.1 Simple Estimate of the Detection Rate 

SN 1987A was discovered optically in the Large Magellanic Cloud, a nearby 
satellite of the Milky Way, at a distance of about 50 kpc (14). Since it was so 
close, it was likely that its neutrinos had interacted in the Kamiokande-II (Kam- 
II) and Irvine-Michigan-Brookhaven (1MB) water-Cerenkov detectors, originally 
intended for proton decay searches. Kam-II reported 12 events and 1MB reported 
8 events, with energies of several tens of MeV, in contemporaneous bursts of about 
10 seconds, occurring a few hours before the optical brightening (15-18). The 
events were largely consistent with being due to inverse beta decay, i>e + P — > 
e+ + n, where the positron is detected, and where its energy is close to that of 
the neutrino. 

Deferring many details, we can simply estimate the (steady) detection rate of 
DSNB neutrinos from the (burst) detection rate of SN 1987A neutrinos, as 
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The detection rate in Kam-II during the burst was ~ 1 s~ . How arc the inputs for 
the dimensionless product of conversion factors related? For SN 1987A, Nsn = 1, 
while there are Nsn ~ 100 neutrino bursts in the universe starting within any 
10-second interval; this accounts for the supernova rate of the Milky Way, 10~^ 
year""*^, the space density of comparable galaxies, ~ 10~^ Mpc~^, and the fact 
that the supernova rate per galaxy was ~ 10 times higher in the past. The 
detector mass M^et of Kam-II was ~ 10 times smaller than that of SK, which is 
22.5 kton (22.5 x 10^ g). The distance D of SN 1987A was 0.050 Mpc, whereas 
for a typical supernova contributing the DSNB, at 2; ~ 1, it is c/Hq ^ 4000 Mpc, 
~ 10^ times farther. Then the time-averaged DSNB detection rate in SK, noting 
changes in Nsn, M^et, and l/D^, respectively, is 

^ ~ (1 s"^) X 100 X 10 X 10^1° ~ lO"'^ s"^ - 3 year-i , (2) 

. dt J DSNB 
which is about the right number. 

The Poisson expectation for the number of neutrinos detected per supernova is 
^ 1, so nearly all supernovas will yield 0, and a rare few unidentified supernovae 
will yield just 1. While the long-term average of the neutrino detection rate 
from Milky Way supernovae is ~ 100 times larger than for the DSNB, it requires 
waiting decades for large bursts from single identified supernovae. Future very 
large detectors may detect more frequent, but much smaller, bursts from identified 
supernovae in the nearest galaxies. Burst detection will be required to study the 
time profile and flavor dependence of supernova neutrino emission. 



3.2 Line of Sight Integral for the DSNB Flux 



The DSNB flux spectrum at Earth is calculated from the neutrino emission per 
supernova ^{En) and the evolving core-collapse supernova rate Rsn{z) as 
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where this form arises as a standard Une-of-sight integral for the radiation inten- 
sity (flux per soHd angle) from a distribution of sources. Neglecting cosmological 
factors, the differential injection rate of neutrinos depends on the rate density (the 
number of neutrinos emitted per volume and time, i.e., the supernova neutrino 
emission times the supernova rate density) times the volume element d^r^dr, 
while the fluxes fall off as l/(47rr^), canceling the r^. The received intensity is 
then just a line integral of the rate density. Since the sources are isotropic on the 
sky, the detection reaction has little directionality, and Earth is transparent to 
neutrinos, the intensity has been integrated over Air to obtain the flux. 

This defines the neutrino number flux spectrum, in units cm~^ MeV~^. To 
obtain the energy flux spectrum, more familiar to photon astronomers, multiply 
through by Ei,. The first term in large brackets is the emission spectrum (p{Ey), 
where a neutrino received at energy E^j was emitted at a higher energy Ei,{l + z)\ 
the prefactor of (1 + z) on the spectrum accounts for the compression of the 
energy scale. The second term is the supernova rate density Rsn{z)- The third 
term is the differential distance, where |dt/dz|~^ = Ho{l + z)[flA + ^Tn{'^ + , 
the cosmological parameters are Hq = 70 km s~^ Mpc~^, = 0.3, S^a = 0.7, 
and their uncertainties can be neglected. 

Above, it is assumed that there is no variation in the supernova neutrino emis- 
sion among different progenitor stars, which follows from the core-collapse cri- 
terion for massive stars and how supernova neutrino mixing effects work, both 
discussed in Section 4.2. It is also assumed that nothing except redshifting hap- 
pens to the neutrinos between the supernovae and Earth, as obscuration is always 
negligible, the neutrino mass eigenstates that propagate do not oscillate, and no 
exotic effects are expected (but such can be tested (19-22)). These are realistic 
assumptions at the precision of the expected counts and energy resolution, and 
can be generalized with appropriate averages. 

3.3 Contemporciry Inputs, Uncertainties, and Detectors 

Throughout, we focus on the effective spectrum at the surface of the super- 
novae, after any neutrino mixing effects in the star, as this is the same observable 
for both burst and DSNB detection. As described in Section 4, the neutrino 
emission in one flavor is characterized by its total energy (time-integrated lumi- 
nosity) and average energy (3.15 T for a Fermi-Dirac spectrum, where T is the 
temperature in units of MeV). These are the parameters to be determined by 
a DSNB measurement, as the supernova neutrino emission is now the most un- 
certain input, and it cannot be measured by astronomical techniques. Nominal 
values are Ep^^tot ~ 5 x 10^^ erg (3 x 10^*^ McV), and Ep^^avg ^ 15 MeV. If a 
black hole forms promptly, the luminosity from the collapsing core prior to being 
abruptly terminated can be higher than for neutron star formation, yielding a 
larger time- integrated emission (24). 

The cosmic supernova rate is now well measured, as shown in Section 5. For 
the supernova rate density, the local value is ~ (1.25±0.25) x 10~^ Mpc~^ year~^, 
and the rate grows by a factor ~ 10 by z = 1. DSNB measurements are also 
sensitive to core collapses that are optically dark but neutrino bright. 

Figure 1 shows the i^g flux spectrum; the SK 2003 limit is < 1.2 cm~^ for 
energies integrated above 19.3 MeV, ruling out spectra above about the middle of 
the T = 8 MeV band. Since the predictions are for generic thermal spectra, they 
could be used for other flavors as well. Equation (3) defines the flux spectrum of 



John F. 




[MeV] 

Figure 1: Predicted DSNB flux spectrum. The spread between different 
labeled temperature bands indicates the uncertainty in the supernova neutrino 
emission; the widths of the bands indicates the uncertainty in the cosmic su- 
pernova rate. Grey shaded regions indicate energy ranges where noted detector 
backgrounds are important. Figure as shown in Ref. (23). 

neutrinos. In Section 6, the event rate spectrum of positrons from the reaction 
I'e + p ^ + n is introduced to allow comparison with the SK data. 

3.4 Theoretical Predictions and Experimental Revolutions 

Theoretical predictions of the DSNB flux go back more than forty years, e.g., 
to brief mentions by Zel'dovich and Guseinov (25) and Ruderman (8), and a 
calculation by Guseinov (26). First examples of detailed calculations, from the 
early 1980s, are Refs. (27-34); unfortunately, the true flux is not as large as in 
these predictions (but see the pessimistic lower limit noted in Ref. (33)). The 
variation between predictions, formerly orders of magnitude, became much less 
beginning in the mid-1990s, in e.g., Refs. (35-45). At roughly the end of that 
era, the DSNB was reviewed in 2004 by Ando and Sato (46). The SK 2003 upper 
limit on the DSNB was a key development, as it directly rules out models with 
fluxes more than a few times larger than present predictions. 

It would be unfair to judge earlier predictions without adjusting their inputs 
to contemporary values. This seems difficult, as a wide variety of different for- 
malisms, inputs, standards of precision, and methods of quoting results were 
used. However, this is easy in the present perspective, first noted in Ref. (47), 
where the principal remaining unknown is the supernova neutrino emission, and 
the assumed values could be easily extracted from earlier papers. In fact, most 
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DSNB predictions did use reasonable supernova neutrino emission parameters 
(for the earliest papers, the existence of i/^ and l'-^ was unknown), even before 
the detection of SN 1987A. There were huge differences in the supernova rate 
density assumed, due to uncertainties that have been dramatically reduced by 
comprehensive star formation and supernova rate measurements over the past 
decade. The results shown in this review cover the range of reasonable allowed 
models for both the supernova neutrino emission and the supernova rate density. 

4 First Ingredient: Supernova Neutrino Emission 

The detection of neutrinos from SN 1987A established that .some supernovae 
produce some neutrinos. Of course, our knowledge is much greater than that - 
there is strong evidence that core-collapse supernovae can be identified optically, 
and that a powerful neutrino burst must accompany core collapse. 

To understand the properties of massive stars and how supernovae explode, it is 
necessary to confirm the ubiquity of neutrino emission from all types of supposed 
core-collapse supernovae and to refine our knowledge of this emission. 

4.1 The Fates of Massive Stars 

Massive stars race towards their doom. The lifetime of a star, mostly the time 
spent burning hydrogen in the core, is ^ 30 (8Mq/M)2-5 X 10^ year, where Mq = 
2.0 X 10^^ g is the solar mass (48). The fusion reactions proceed in stages that 
depend on temperature, overcoming the Coulomb barriers for successively heavier 
elements only once the preceding elements are exhausted and the core contracts. 
Only stars with mass greater than about 8Mq will become hot and dense enough 
in the core to burn past carbon and oxygen to iron (49,50). A Chandrasekhar 
instability sets in once about 1.4Mq of iron has been created, since no further 
nuclear reactions can generate energy and electron degeneracy pressure provides 
insufficient support. Photo-nuclear reactions then endothermically destroy iron, 
leading to further collapse, which continues until the core reaches miclcar densities 
and cannot be compressed further. The subsequent bounce leads to the formation 
of an outgoing shock, which removes the envelope of the star and causes the 
optical supernova, which is bright for several months. What will be left behind is 
a neutron star or black hole from the core and a shell remnant from the envelope. 

Astronomical characterizations of supernovae are based on their optical prop- 
erties (51). Supernovae must be distinguished by their spectra, as their light 
curves alone are not usually enough. Many lines of evidence indicate that Types 
II, lb, and Ic supernovae arise from the core collapse of massive stars, differing 
in the amount of progenitor envelope mass loss before explosion. Type II su- 
pernovae are several times more common than Types lb and Ic, and the latter 
two seem to arise from ranges of successively larger progenitor mass. All should 
have comparable emission of all flavors of neutrinos and antineutrinos, due to the 
mechanism of core collapse. Type la supernovae, which are also several times 
less common than Type II supernovae, seem to be powered by the runaway ther- 
monuclear burning of a carbon and oxygen white dwarf in a binary system into 
iron; only relatively modest emission of low-energy electron neutrinos is expected. 
Detecting neutrinos will test the mapping between optical supernova types and 
underlying physical mechanisms (52). 
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4.2 Core-Collapse Supernova Neutrino Emission 

The basic picture of core collapse is well known (53). When the core collapses, 
there is a huge change in its gravitational potential energy. For a constant density 
profile, the change in this self-energy is 



where the number follows from the mass and radius of a typical neutron star, 
IAMq and 10 km (the subtracted term is negligible since the initial core radius is 
thousands of kilometers), and this is an incredible 10% of Mc^ for the core (54- 
56). The mass of the collapsed core varies little with the mass of the progenitor 
star (57, 58) . Since the proto-neutron star is near nuclear density, no particles 
except neutrinos can escape, and core collapse is inevitably accompanied by a 
huge burst of neutrinos, carrying away nearly the whole change in gravitational 
potential energy. The neutrinos, which begin leaving the core at the instant of 
collapse, arrive at Earth before the light, which is delayed by at least the time 
it takes the shock to propagate through the envelope, typically hours. This was 
observed for SN 1987A, and is a crucial confirmation of the basic picture. 

It is thought that each of the six flavors of neutrinos and antineutrinos carries 
away a comparable fraction of the total energy. In numerical simulations of 
supernovae, this equipartition is only approximate, and the sense in which it 
might be violated in nature is not known (59-62) . Most of the neutrino emission 
is in pairs, i/e + Vg, I'n + i'ix, and Ut + i^t, that are produced by weak interactions in 
the hot and dense matter in the core. The reaction e~ + p ^ v^ + n converts the 
nearly equal numbers of neutrons and protons from nuclei into the neutron-rich 
matter of the neutron star, producing an excess Ve flux, but this is a ~ 10% effect. 

The proto-neutron star density is so high that even neutrinos do not read- 
ily escape, and instead diffuse out over several seconds, and are thus emitted 
with a quasi-thermal spectrum characteristic of the surface of last scattering. 
If the proto-neutron star were semi-transparent, the neutrinos would leave on a 
timcscale of ~ 1 millisecond, with energies of 100 MeV. Taking opacity into ac- 
count, and assuming blackbody emission, L ^ Etot/T ~ B?T'^, the increase in the 
emission time by ~ 10^ corresponds to the decrease in neutrino average energy 
by ~ 10. A long emission timescalc and moderate temperature were observed for 
SN 1987A, confirming the high density of the proto-neutron star. 

The opacities of neutrino fiavors with matter vary, leading to different neutrino 
temperatures, since the matter temperature of the proto-neutron star falls with 
radius. The interactions of u^^ v^-^ Vr arc only through the weak neutral 
current, since their energies are too low to produce the corresponding charged 
leptons. For and Ve-, there are also interactions through the weak charged 
current, so the cross sections are larger; for i/g, the coupling constants are larger 
and there are more neutron than proton targets. Oft-quoted estimates for the 
temperatures are T = 8 MeV for v^, v-r, I'lxi T = 5 MeV for Pg, and T = 4 
MeV for v^. More recent theoretical work indicates that the temperatures, and 
the splittings between them, arc less (64). 

Neutrino mixing effects relate the initial neutrino spectra at the opaque emit- 
ting surface of the proto-neutron star to those outside the supernova. While the 
opacity is low outside the neutron star, the high density of matter and even other 
neutrinos significantly affects neutrino mixing (65-68). The passage through high 
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Figure 2: Time-integrated effective v'e spectrum from SN 1987A. The 

curves indicate results reconstructed using only the Kam-II data, the 1MB data, 
or their combination. The grey shaded region is a thermal spectrum with param- 
eters Ep^^tot ~ 5 X 10^^ erg and Ep^^avg ~ 15 MeV. Figure as shown in Ref. (63). 

densities leaves the neutrinos in incoherent mass eigenstates that do not oscillate 
between the supernova and Earth. Mixing effects in Earth are suppressed by the 
isotropy of the sources and the detection cross section (69). 

Neither the expected initial spectra nor the neutrino mixing effects are un- 
derstood well enough, so we use the observable effective spectra. For u^, the 
effective time-integrated spectrum after mixing, assumed to approximately be of 
the Fermi-Dirac form, in units MeV~^, is 



where the total energy is Ep^^tot and the average energy is Ep^^avg = 3.15 T. These 
parameters are not fixed from theory, but will be determined from experiment, 
and will test supernova simulations and neutrino mixing scenarios. 

4.3 Comparison to Supernova 1987A 

While the SN 1987A data were sparse, and the detectors were not optimized for 
that purpose, these data are the most direct evidence we have about supernova 
neutrino emission. A reconstruction of the effective i>e spectrum from the Kam-II 
and 1MB data is shown in Figure 2; it is in reasonable agreement with expecta- 
tions. Higher temperatures are favored by theoretical models including neutrino 
mixing, as well as by observations of nucleosynthesis yields that depend on high- 





John F. Beacom 



11 



energy neutrinos (70-72). It remains to be seen if SN 1987A was different from 
an average supernova, and more data are needed, especially at high energies. 

The SN 1987A data are an important input for estimating the DSNB spec- 
trum (40,45), though care must be taken to note that the high-energy SN 1987A 
neutrinos are the most relevant for predicting the DSNB (63). The SK 2003 
search used a positron energy threshold of 18 MeV, and future searches may 
reach 10 MeV; redshift effects enhance the importance of high-energy emission. 

5 Second Ingredient: Cosmic Supernova Rate 

Supernovac arc infrequent in the Milky Way, but not in the universe. Since 
massive stars have lifetimes that are very short on cosmological timescales, their 
cosmic birth and death rates are exactly equal. The cosmic star formation rate 
has been measured precisely, using a great variety of techniques, principally based 
on the emission of massive stars. This is enough to accurately give the rate of core 
collapses, and it is supplemented by direct measurements of the optical supernova 
rate, which have lower precision, but are in good agreement. 

5.1 Measured Star Formation Rate 

Star formation rate measurements do not literally probe star formation, but begin 
with the total luminosity of the stars in a galaxy. If these stars were all of the 
same known mass, so that their luminosity and lifetime were given by stellar 
evolution theory, then the star formation rate would be determined by the total 
mass of these stars divided by their lifetime; this is the birth (and death) rate 
needed to keep the number of stars in equilibrium. The main star-formation 
rate indicators measure the emission from massive stars that are short-lived on 
galactic timescales, so the equilibrium assumption should be good on average. 

A wide range of stellar masses are always present, complicating this proce- 
dure, but the degeneracies can be broken, since stars of increasing mass have 
higher temperatures and much higher luminosities. Massive stars can be iso- 
lated by measurements of their continuum flux or of nebular emission lines 
from atomic recombinations in the surrounding gas they ionize. The distri- 
bution of stellar masses formed when star formation occurs is assumed to fol- 
low a universal initial mass function (IMF), e.g., the conventional Salpeter IMF 
scales as 7/>(M) = dn/dM oc M~^-^^ for stellar masses M between O.IMq and 
lOOM© (73,74); more realistic IMFs modify the low-mass end. With an IMF, 
calculated stellar and nebular emission spectra, and careful consideration of the 
effects of obscuration and re-radiation, the star formation rate can be determined 
from a galaxy luminosity spectrum (75). The high-mass star formation rate is 
measured, but the formation rate of low-mass stars has to be deduced assuming 
an IMF. This introduces a factor-two range in the normalization of the total star 
formation rate, due to the low-mass IMF uncertainties. 

The comoving star formation rate density, in units Mq Mpc~^ year~^, is about 
ten times larger at z = 1 than it is today, at z = 0. Comoving coordinates are 
those for which the expansion of the universe has been removed; crudely, this 
defines the star formation rate per average galaxy. Aside from the normalization 
uncertainty due to choice of IMF, the precision of the combined data is at the 
tens-of-percent level for 2; < 1 (23, 76). From 2; = 1 to at least 2^4 — 5, the 
star formation rate is nearly flat, with a precision becoming as poor as a factor 
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of two, until a more uncertain decline at higher redshifts, e.g., Refs. (77-82). 

These results, especially those in the z < 1-2 range relevant for the DSNB, 
are quite robust, since they have been measured by many different groups using 
different instruments and techniques, and extensive efforts have been made to 
calibrate the observables and their corrections for dust obscuration and incom- 
pleteness. The star formation rate data are also in good agreement with mea- 
surements of the cxtragalactic background light, which, analogous to the DSNB, 
is an integral measure of stellar emission over redshifts (23). These and other 
data favor an IMF of somewhat shallower slope than the Salpeter IMF. 



5.2 Predicted Supernova Rate 

To predict the core-collapse supernova rate density, Rsn{z), in units Mpc~^ 
year~^, from the star formation rate density, Rsf{z), we use 

Rsn{z) = Rsf{z) —TKFi — , (6) 

^ ' ^ ' ffM'^{M)dM 143M0 ^ ' 

where the upper integral gives the number of stars that lead to core collapse and 
the lower integral gives the total mass in stars, and a Salpeter IMF has been used. 
Importantly, for Rsn{z) the above-mentioned normalization uncertainty due to 
the choice of IMF goes away, because nearly the same massive stars source the 
star formation rate measurements and lead to core-collapse supernovae. For an 
IMF with a shallower slope, Rsf{z) is smaller, but the ratio of integrals is larger, 
canceling in Rsn{z) to the few-percent level (23). 

In the upper integral, the choice of upper limit of integration is unimportant, 
as long as it is large; then the integral scales as {SAlQ/MiowerY'^^ for a Salpeter 
IMF. Theoretical and indirect empirical work indicates that the minimum mass 
star that leads to core collapse is about 8Mq. This is now supported by direct 
observations. Using archival imaging from the Hubble Space Telescope and other 
sources, it has been possible to identify the stellar progenitors of some supernovae 
in nearby galaxies, e.g., Refs. (83-86). The luminosity of the star determines its 
mass, and a fit to progenitor data gives a minimum mass to lead to a core-collapse 
supernova of (8.5 it 1.5) Mq (87). In addition, the largest estimate so far of the 
maximum mass star that ends instead as a white dwarf is ~ 7Mq (88,89). 

There is therefore a precise prediction of the cosmic core-collapse rate, based 
on strong evidence from star formation rate data and direct confirmation of the 
minimum stellar mass to lead to a core-collapse supernova. This is based on 
observations of living massive stars, and does not differentiate whether they will 
die as visible supernovae or optically-dark collapses. Since those outcomes have 
comparable neutrino emission, this data is enough to precisely predict the DSNB. 



5.3 Measured Supernova Rate 

As noted, the relationship between star formation and supernova rates is just 
a simple conversion with small uncertainties on the result. Since the redshift 
evolution of the star formation rate is well measured, the evolution of the super- 
nova rate is then precisely defined, and the conversion predicted above could be 
checked at a single redshift. 

Measuring supernova rate densities is not easy, since supernovae are infrequent 
in individual galaxies. Core-collapse supernovae (Types II, lb, Ic) are fainter and 
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Figure 3: Predicted and measured cosmic core-collapse supernova rates. 

In the upper panel, the band and its width indicate the absolute prediction and its 
uncertainty from the star formation rate data alone (shown, scaled down, with 
faint grey points). The most reliable supernova rate measurements are shown 
with solid symbols; those with open symbols must be less complete. In the lower 
panel, rate ratios are shown, which supports this. Figure as shown in Ref. (23). 



less uniform in their optical properties than the thermonuclear supernovae (Type 
la) used as distance indicators, and have received much less attention, despite 
being more common. 

At high redshift, repeat imaging of a volume can find supernovae largely inde- 
pendent of their host galaxies. This directly measures the quantity of interest. At 
low redshift, galaxies from a catalog can be monitored. Since this finds the galax- 
ies first, and the supernovae second, it is subject to incompleteness due to faint 
galaxies, and one must carefully define the global volume and time monitored 
from those for each galaxy. While these corrections are made, catalog surveys 
tend to find lower supernova rates and less extreme supernovae than volume sur- 
veys. Lastly, supernova rates estimated using supernovae that were not collected 
in a systematic way should be taken only as lower bounds. 

The predicted supernova rate density is shown in Figure 3. There is very good 
agreement with the most reliable data, despite the conversion being predicted 
absolutely, and not fit. The uncertainty on the star formation rate is taken 
as 20% on Rsf{z = 0). The predicted supernova rate and its uncertainty are 
Rsn{z = 0) = (1.25 ± 0.25) X 10"^ Mpc^^ year"! (recall Section 3.1). At 
redshift zero, the lower limits on the rates are already quite high, matching the 
prediction well; the 10 Mpc was defined very conservatively, and this and the 
28 Mpc point are known to be incomplete (52,87). Therefore, supernova rate 
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measurements lower than this band must be less complete, since the shape is 
fixed from the star formation rate. At high redshift, the updated Dahlen et al. 
points are in excellent agreement; these were measured volumetrically (90,91), 
and are larger than points that were not (92-94). [After Ref. (23) was published, 
the Dahlen et al. points decreased slightly, but are still consistent with the band 
shown in Figure 3 (T. Dahlen, private communication, 2010).] That faint core- 
collapse supernovac were missed is supported by data on their rate relative to the 
brighter Type la supernovae. Note that the discovery rates of supernovae have 
dramatically increased with time, with advances in technology and interest. 

The agreement of the predicted and measured supernova rates strongly sup- 
ports this input for the DSNB calculation. While a significant fraction, 50%, 
of dark collapses is allowed, this is not required (~ 10% is expected). These col- 
lapses might be intrinsically faint, truly dark, or just obscured. In any case, their 
existence can only raise the level of the DSNB. In the near future, supernova rate 
measurements will greatly improve, and their comparison with star formation rate 
data will allow more precise tests of the fraction of dark collapses (23,95-97). 

6 Third Ingredient: Neutrino Interactions and Detectors 

The detection probability depends on the product of the DSNB flux and the 
effective area of the detector. Since SK (and Earth) is nearly transparent to 
neutrinos, the effective area is the neutrino-proton cross section times the number 
of target protons, independent of the shape or orientation of the detector. 

While supernovae make all flavors of neutrinos and antineutrinos, the best 
prospects are for i>e in SK, detected hy i>f, + p ^ + n, as it has both large size 
and low background rates. The properties of neutrinos are now measured well, 
and exotic scenarios that would affect DSNB detection, e.g., large mixing with 
sterile species, are no longer favored. 

6.1 Detection of DSNB Electron Antineutrinos 

Electron antineutrinos are detected by the inverse beta decay reaction on free 
protons, i>e + p ^ e^ + n. Here "free" protons means hydrogen nuclei, and not 
the protons bound in heavier nuclei, for which nuclear binding effects suppress 
interactions at low energies; the atomic ionization state is irrelevant. 
The total cross section for inverse beta decay is 

a{E^) = [0.0952 x 10"^^ ^^^2 _ ^ ^ MeV)^] (1 - lE^/Mp) , (7) 

for neutrino energies above threshold, Ey > 1.8 MeV, and where Mp is the proton 
mass (for c = 1). In the high-energy range, the rising cross section helps offset the 
falling flux spectrum. This expression, which is good to the several-percent level 
in the relevant energy range, neglects the small positron mass but does include 
the few-tens-percent recoil-order corrections, shown outside the square brackets. 
At the same orders, the kinematic relationship is 

Ee = [{E, - 1.3 MeV)] (1 - EjMp) , (8) 

where the isotropy of the neutrino angular distribution has been taken into ac- 
count. The recoil-order correction to the kinematics can be neglected relative to 
the detector energy resolution. Full expressions for the cross section and kine- 
matics are given in Ref. (98,99). 



John F. Beacom 



15 



6.2 The Super-Kamiokande Detector 

The fiducial volume of SK contains 22.5 kton of water, and hence 1.5 x 10^^ free 
protons (100). To reduce backgrounds, SK is under more than 1 km of rock, is 
heavily shielded by separate outer and inner buffer regions, and was constructed 
to maintain extreme purity. SK is large enough to have a few DSNB interactions 
per year, and already has low background rates. These facts are fortuitous, as 
the SK design was optimized not for the DSNB, but for measurements of proton 
decay, and atmospheric, accelerator, solar, and Milky Way supernova neutrinos. 

In SK, relativistic charged particles are detected by the cones of optical Cerenkov 
light produced as they lose energy; the Cerenkov process is a small component of 
the energy loss rate, but is the only one that SK detects directly. The Cerenkov 
light travels through water to the photomultiplier tubes on the walls, which view 
a homogeneous volume of transparent water. From the patterns of received pho- 
tons, the position, direction, energy, and identity of charged particles are deter- 
mined. SK has excellent detection efficiency for the energies considered here. 



7 Framework Redux: DSNB Detection Spectrum at Earth 



Now that we have reviewed the three ingredients - the supernova neutrino emis- 
sion, cosmic supernova rate, and detector capabilities - we return to the frame- 
work for the DSNB detection spectrum at Earth. 

The event rate spectrum, in units s~^ MeV"-*^, where the detection cross section 
(j{Ey), the number of proton targets Np, and the shift between neutrino energy 
Ey and positron energy — — 1.3 MeV have been accounted for, is 



f-OD 

(Ee) = Npa{E,) {l + z)ip[E,il + z)] Rsn{z) 



cdt 



dz 



dz 



■ (9) 



For relevant energies, those above 10 MeV, the effective upper limit on the line of 
sight integral is ^; ~ 1-2, as beyond there the falling dt/dz is no longer compen- 
sated by a rising Rsn{z) and the larger emission energies required for the same 
observed energy probe further into the high-energy tail of the emission spectrum. 

In Figure 4, we show predictions for the DSNB detection spectrum, taking 
up-to-date inputs and their uncertainties into account. The range due to differ- 
ent emission parameters, especially at high energies, indicates the uncertainties 
to be reduced by measuring the DSNB. The uncertainty on the supernova rate 
is already modest and will soon be smaller. The detection reaction is inverse 
beta decay in SK, which is well understood. This predicted event rate spectrum 
can be directly compared with measured data, which are a combination of DSNB 
signal and detector backgrounds. The DSNB is especially sensitive to the im- 
portant high-energy part of the emission spectrum. The predictions for different 
temperatures converge near 10 MeV, so a measurement there would probe Ej^^^tot 
alone; the falloff of the spectrum would then probe Ep^^avg alone. 



8 Super-Kamiokande 2003 Upper Limit on the DSNB 

It takes extraordinary efforts to operate a gigantic detector like SK with low 
background rates at low energies. The SK DSNB limit (101) improved upon the 
Kam-II limit (102) by a factor ~ 100, establishing the only experimental limit 
within range of contemporary predictions, an extremely important achievement. 
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Figure 4: Predicted DSNB z^g event rate spectrum in positron energy. 

The labeled bands and their widths are as in Figure 1. Integrated event rates are 
tabulated in Ref. (23). The 2003 exclusion from SK is shown; it is largely inde- 
pendent of the assumed temperature. The energy range of the irreducible reactor 
background is shaded; backgrounds at higher energies depend on if gadolinium is 
added to SK. Figure as shown in Ref. (23). 

8.1 Detector Backgrounds and DSNB Flux Limit 

For the SK 2003 DSNB search, only visible positron or electron (they cannot 
be distinguished) energies above 18 MeV were used, where detector background 
rates, due to atmospheric neutrinos, are relatively low. Nothing beyond expec- 
tations was found in this pioneering search, including any surprise problems that 
would impede this or improved future searches. 

To understand the atmospheric neutrino backgrounds, it is essential to use the 
event rate spectrum in visible energy. The flux spectrum alone is misleading, as 
here the visible energy is quite different than the neutrino energy, unlike for the 
DSNB signal. The backgrounds, shown in Figure 5 with visible energies of 18- 
82 MeV, arise from neutrinos with energies up to ~ 250 MeV that invisibly enter 
the detector and have charged-current interactions, mostly with oxygen nuclei, 
inside the fiducial volume. 

Atmospheric z/^ and can produce nonrelativistic /i~ and /x"*', which do not 
produce Cerenkov light as they lose energy. When these decay at rest, the elec- 
trons and positrons produced are relativistic, and their spectrum is the bump in 
Figure 5. This is the dominant background, but its shape is fixed and its nor- 
malization can be measured above the energy range of the DSNB. Atmospheric 
Vf, and Ve can produce e~ and e"*" with energies well below the neutrino energy, 
due to nuclear effects. These produce the rise at high energies in Figure 5. 
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Figure 5: Results of the SK 2003 DSNB z/g search. The efficiency-corrected 
measured data are shown with points and error bars. The expected total atmo- 
spheric neutrino background is shown by the thick solid line, and its components 
by the thin solid lines. The largest allowed DSNB signal is shown by the shaded 
region added to the atmospheric background. Figure adapted from Ref. (101): 
here dN/dEg is shown, to be directly compared in 1-MeV steps to Figure 4; to 
recover the dN results shown in Ref. (101), multiply the bin values by 4 MeV. 



A DSNB signal would have appeared as an excess at low energies (compare to 
Figure 4); none was seen, and a limit of (p{Ei^ > 19.3 MeV) < 1.2 cm~^ at 
the 90% confidence level was set (101). 

8.2 Calculated Limits on Supernova Neutrino Emission 

A DSNB limit expressed as a flux integrated above a certain energy, while a good 
starting point, is model-dependent and makes comparisons difficult. One variable 
is not enough to characterize DSNB spectra, but the two natural variables of total 
energy Ep^^tot and average energy Ep^^avg for the Ue flavor are sufficient (47). These 
directly control our ignorance of the effective supernova neutrino emission, since 
the cosmic supernova rate is now known well enough. This reinterpretation of 
the SK DSNB limit, showing how the two emission parameters can compensate 
each other, is shown in Figure 6. Any change in the assumed normalization 
Rsn{z = 0) corresponds to a simple change in Ep^^tot alone, and uncertainties on 
the evolution Rsn{z) / Rsn{z = 0) are small. 

Excitingly, the sensitivity is near the parameters of supernova models (59-62) 
and the uncertain fits to the SN 1987A data (103). A two-parameter limit on the 
effective Ve emission is simple for experimentalists to implement and corresponds 
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Figure 6: Constraints on the effective supernova emission parameters. 

The grey excluded region follows from the null result of the SK 2003 DSNB search. 
Approximate regions from fits to the SN 1987A data are shown, e.g., Ref. (103). 
The point D indicates the nominal parameter values. Figure as shown in Ref. (47). 

to the minimal theory. If it is proven that there is large variety in supernova neu- 
trino emission, large neutrino mixing effects, or exotic effects, then appropriate 
theoretical averages should be compared to this direct experimental limit. 

9 Pathways to Discovery 

Progress on experimental sensitivity will have a decisive impact, as the expected 
signal is close to current limits. Astronomical observations and theoretical work 
will refine and enhance these prospects. 

There are two striking features of the SK 2003 results. First, in the energy 
range used, the search is background-limited, so the sensitivity will not improve 
linearly with exposure time, but naively only with the square root. Second, going 
to lower energies would give much more DSNB signal and much less atmospheric 
neutrino background; however, there is a large background rate due to the beta 
decays of nuclei produced by spallation interactions of cosmic-ray muons. To 
detect the DSNB, SK therefore needs to reduce backgrounds and extend the 
search to lower energies. But how? 

9.1 Proposal for Background Reduction in SK 

Since the DSNB detection reaction is z^e + P ^ e"^ -|- n, while the detector back- 
ground reactions mostly do not produce neutrons, it is obvious that neutron 
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Figure 7: Expected detection rates in SK with dissolved gadolinium. 

The DSNB signal could be cleanly detected at a reasonable rate (see Figure 4 
for updated DSNB predictions and uncertainties). Figure as shown in Ref. (104). 



detection is the key to separating signal from backgrounds. The neutron loses 
energy by elastic scattering and, at present, captures on a free proton, producing 
a 2.2-MeV gamma ray that leads to an insufficient detection signal in SK. 

The idea of positron-neutron coincidence detection goes back to Reines and 
Cowan, and it is well known how to implement this in oil-and-scintillator de- 
tectors, which have much higher light yields than water- Cerenkov detectors. For 
example, with the suspension of a small concentration of a gadolinium compound 
in the oil, neutrons will capture on gadolinium, due to its enormous cross section, 
producing an easily detectable 8- MeV gamma-ray signal. The neutron is detected 
a long tens of microseconds after the positron, at nearly the same position. 

However, until Ref. (104), it had not been shown that neutron detection with 
dissolved gadolinium in a water-Cerenkov detector could be feasible. The key 
points are that soluble compounds exist, e.g., GdCls and Gd2(S04)3, and that 
the gamma rays would lead to a detectable signal in SK, due to its high photo- 
multiplier coverage. That these and many other technical concerns had positive 
answers was quite surprising. The ability to detect neutrons would give unprece- 
dented capabilities to separate signals and backgrounds in very large detectors. 
As described in Ref. (104), this would reduce the atmospheric neutrino back- 
grounds and, more importantly, would dramatically reduce the large spallation 
radioactivity backgrounds, allowing SK to reduce its DSNB analysis threshold to 
about 10 MeV, below which the reactor background is overwhelming. 

Figure 7 shows that there would be excellent prospects for the DSNB if SK had 
a ~ 0.1% concentration of dissolved gadolinium. There is no known astrophysical 
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source that can produce a comparable high-energy Aux, leaving a detection 
window for the DSNB. If an exotic process converts a small fraction of solar 
to then the discovery of this new physics could partially get in the way of the 
DSNB, but we should have such problems (104, 105). 

9.2 Research and Development for Gadolinium in SK 

Prom the time of Ref. (104), Vagins, aided by other SK experimentalists, has 
been intensively researching the practicalities of introducing a dissolved gadolin- 
ium compound into SK, e.g., Refs. (106-109). Our initial investigations, and the 
more extensive research that followed, have uncovered no insurmountable prob- 
lems regarding the infrastructure required to manage the gadolinium, the effects 
on the detector and other physics searches, or the positron-neutron coincidence 
detection. However, it has been determined that it is not advisable to pack a 
kilogram of white GdCla powder in carry-on luggage for an international flight. 

The progress on the comprehensive research and development program, and 
the seriousness with which the SK Collaboration is considering adding gadolin- 
ium, can be illustrated with two recent examples. First, a small container of 
gadolinium-loaded water and a neutron source were lowered into SK, and it was 
found that the detectable signal of neutron capture on gadolinium works as ex- 
pected (110) (see also Ref. (111)). Second, a new 200-ton "miniature" (~ 1% 
scale) of SK with dissolved gadolinium is being built underground near SK to 
test all aspects of a gadolinium-loaded water-Cerenkov detector (112). 

9.3 Prospects for Other Detectors and Flavors 

Oil-and-scintillator detectors the size of SK would have reduced atmospheric 
backgrounds, but the same DSNB signal rate and a reactor background that is 
relatively large at any location on Earth (113,114). Ultimately, a high-statistics 
measurement of DSNB will be needed to fully develop the science possibilities, 
probably with a water-Cerenkov detector at least ten times larger than SK, as is 
being considered for other purposes (115,116); the prospects for a high-statistics 
DSNB measurement have been explored (117-121). 

Direct DSNB limits on other neutrino flavors are weaker than for i?e in SK. With 
all their data analyzed, the Sudbury Neutrino Observatory sensitivity to DSNB 
fe will only be several times worse (122,123), comparable to indirect limits (124), 
but limits for other flavors are worse by orders of magnitude (125, 126). Large 
liquid-argon detectors are promising for DSNB Ug (127,128). 

10 Conclusions 

Neutrino astronomy can revolutionize astrophysics and physics, due to its unique 
abilities to see unseen processes and to probe subtle new particle interactions. 
Stellar fusion reactions and the core-collapse mechanism were revealed by the 
detection of neutrinos from the Sun and the nearby SN 1987A. Solar neutrinos 
were crucial to the discovery of neutrino mixing, and supernova neutrinos probed 
exotic physics far beyond the reach of laboratory experiments. 

However, the great promise of neutrino astronomy remains largely unfulfilled, 
as no other neutrino sources have yet been detected. Recall the courage of the 
pioneers of this field, who, more than forty years ago, had neither adequate 
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detectors nor accurate astrophysics to fortify their hopes. We have advantages of 
many orders of magnitude on each, and can see the end of the desert. To make 
progress on our understanding of supernovae and their role in astrophysics, as 
well as of the neutrinos themselves, we must detect more supernova neutrinos. 

Fortunately, first detection of the DSNB is within reach. Why are we so sure? 

• The framework for calculating the DSNB is solid. 

It is built on simple principles of cosmology, well tested with photons. 

• Supernova neutrino emission has been measured from SN 1987A. 

Theory and indirect evidence suggest it should typically be larger. 

• The uncertainty on the cosmic supernova rate is modest. 

Observations of massive stars and supernovae arc precise and qiiickly improving. 

• All aspects of DSNB signal detection are understood. 

Neutrino interactions and detector properties have each been well measured. 

Therefore, even including uncertainties, there is a robust prediction that there 

must he a few DSNB neutrino events per year in SK. These events are hidden by 
backgrounds - what about the experimental prospects for revealing them? 

• SK has already shown sensitivity quite close to these predictions. 

Importantly, SK has controlled detector backgrounds nearly to the required level. 

• There are good prospects for improving the SK DSNB sensitivity. 

With gadolinium, SK would reduce its backgrounds and energy threshold. 

Therefore, the DSNB could be discovered in SK. With sufficient detector back- 
ground redTiction, even a few signal events will be statistically significant and 
scientifically important. If the DSNB is not found, there must be surprising new 
physics, requiring an upheaval in astrophysics or exotic new particle physics. 

The enduring value of a DSNB measurement will come from its comparison to 
other observations and theory, which include the SN 1987A data and numerical 
simulations of core collapse. The SK 2003 limit is already an important new 
result, as it excludes neutrino emission as large as often assumed, and is only a 
factor 2-4 away from more conservative cases. Significant developments in the 
theoretical predictions are also eagerly anticipated, as, after more than forty years 
of work, recent results have increasingly realistic physics and accurate calcula- 
tions, show explosions in some cases, and some are able to integrate to several 
seconds. These models will provide a rich context for interpreting DSNB results 
and, eventually, to novel tests of neutrino properties and other new physics. 

In the next several years, we will gain a newly deep and comprehensive under- 
standing of supernovae and their neutrinos, following from the above improve- 
ments converging with unprecedented electromagnetic data, results from gravi- 
tational wave observatories, and, with a little luck, a Milky Way neutrino burst. 
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